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Using the non- equilibrium Green's function (NEGF) formalism within the sequential regime, we 
studied ultrahigh spin thermopower and pure spin current in single-molecule magnet(SMM), which is 
attached to nonmagnetic metal wires with spin bias and angle {0) between the easy axis of SMM and the spin 
orientation in the electrodes. A pure spin current can be generated by tuning the gate voltage and 

1 3 

temperature difference with finite spin bias and the arbitrary angle except of 9 = - n, - n. In the linear 

regime, large thermopower can be obtained by modifying V g and the angles {0). These results are useful in 
fabricating and advantaging SMM devices based on spin caloritronics. 



Studies on nanoscale thermoelectric devices have attracted much attention during the past a few years 1-4 . 
It is well accepted that nanoscale materials may provide an opening for the thermoelectricity in 
meeting the challenge of being a sustainable energy source 5 . Huge deviation from the Wiedemann- 
Franz law 6,7 in the nanostructure materials 5 makes new opportunities for investigating novel thermoelectric 
devices with high efficiency 8,9 . Specially, spin caloritronics (spin Seebeck effect) was observed by Uchida et 
al 10,11 . They found that the spin -polarized currents (If and l\) can be induced by a temperature gradient and 
flow in opposite directions. These wonderful discoveries strongly promote research on new energy of 
thermoelectricity 1213 . 

A single-molecule magnet (SMM) is a typical nanoscale material. In experiments, controlling the molecu- 
lar spin 14 and measuring thermopowers of molecule 15 " 20 have been realised by directly using a scanning 
tunneling microscope. The spin-dependent transport properties, such as tunneling magnetoresistance(TMR) 
and spin Seebeck effect, were investigated in the sequential, cotunneling, and Kondo regimes using Wilson's 
numerical renormalization group and quantum master equation 21 " 25 . Many fantastic phenomena have been 
found in the experimental and theoretical studies, including negative differential conductance 26 ' 27 , Berry 
phase blockade 28 , the magnetization of SMM controlled by spin-bias and thermal spin-transfer torque 29 . A 
SMM in a single spin state is necessary for generating the pure spin current 2 1,22,29 without the magnetic field 
or magnetic electrodes. Meanwhile, it implies that the system temperature is limited by the blocking 
temperature of SMM (T B ). When the symmetry of spin in the leads is broken, the angle (9) between 
the easy axis of SMM and the spin orientation in the electrodes will influence the transport properties in 
the SMM devices. Specially, spin-bias 30,31 and this angle (9) are important and crucial on thermoelectric 
effect. 

In this paper, we theoretically investigate the thermoelectric effects of a sandwich structure of NM/SMM/ 
NM with spin-bias 29,32,33 and angles (9) between the easy axis of SMM and the spin orientation in the 
electrodes. We show that, in this system, pure spin currents are observed even though the system temper- 
ature is higher than the blocking temperature due to the spin symmetry broken by spin bias. In the linear 
regime, both thermopower and figure of merit are dependent on the angle and spin bias. It's worth noting 
that the angle plays a critical role on generating spin thermopowers. The figure of merit could tend to infinity 
by tuning the voltage gate at special angles, which implies that this system has an ultrahigh thermoelectric 
efficiency. 
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Results 

Effective hamiltonian. The general Hamiltonian is expressed as 2 ' 
H = Hi ea ds + H SM M + H„ in which 

ueL,R,a'e{ + ,-},k 

H S MM=(to-eV g ) ^2 n a + Un^ni-K 2 S 2 z -Js'S 



0<x t . t 

— C i i.— sin — 
2 a ' + ' fc 2 



^0£ f ^0£ f 



+ H.C 



Hieads describes the free electrons in two leads, with c£ G , k (c^ 0 >k) being 
the creation (annihilation) operator for a continuous state in the 
oceL, R lead with the energy ^ a ^ k and spin index a = + ( — ), 
which denotes spin-majority (spin-minority) electrons. In this 
paper, wideband approximation is adopted and the density of 
states of the leads does not depend on the energy of the two leads. 
The chemical potential of a lead is defined as ii% = +(eV 0l /2 + 



S a >eV sa /2) with eV l 

,L(R) 



u L } R) + u L W 



-AC 



-, and ^ = + l(-l) for </=+(-). V=V L -V R is 

the voltage and V s = V sL — V sR is the spin voltage. P a denotes the 
polarization of a lead and is defined as P a = (p a+ — p a _)/ 
(p a+ + p a _). H^mm denotes the molecular degrees of freedom, in 
which n G = d\d G and d\(d G ) is the creation (annihilation) operators 
for the LUMO. e 0 is the single-electron energy of the LUMO level, 
which is tuned by a gate voltage V g . U is the on-site Coulomb 
repulsion. / describes the Hund's rule coupling between the giant 
spin S of SMM and the electron spin in the LUMO, and parameter 
K 2 is the easy-axis anisotropy of SMM. H t describes the tunneling 
between the LUMO of SMM and the electrodes, and 9 a denotes the 
angle between the spin orientation of lead- a and the easy- axis of the 
SMM (as z-axis). 

In the following, we turn to numerical calculations with para- 
meters: S = 2, / = 0.1 meV y K 2 = 0.04 meV y U = 1.0 meV and 
e 0 = 0.2 meV. The tunneling parameters are set to r = T L = T R = 
0.001 meV. The properties of the leads are set to P R = P L = 0. 
Conventionally, I c and I s are defined as charge current and spin 
current respectively, and we set e = k B = fr = 1 and the thermopower 

and current are scaled in the unit of I 0 = - T. We can find all of the 

thermopowers are symmetric about 9 = n because of the spatial 
symmetry of the sandwich structure. 

Transport properties. First, we consider that the left electrode is 
nonmagnetic with V s = 0.01 meV and V = 0 meV. The system 
temperature is lower than the anisotropy-induced energy barrier 
(K 2 S 2 = 0A6meV). Fi gure la and Figure lb show I c and I s as 
function of 9 for different values of V g with AT = 0.0002 meV, 
respectively. In this case, we can find that I s is almost ten times of 
I c . The maximum or minimum value of I c and I s depends on the gate 
voltage V g and 9, but the positions of these extremums only depend 
on the Vg. When 9 = n/2,3n/2, I s is exactly equal to zero due to the 
coefficient A L aa , = 1/2, which leads 7| = . Moreover, we show I c and 
I s as a function of V g with two types of AT at 9 = 0 in Figure le. One 
can find that I c is extremely sensitive to the temperature difference. 
However, I s only has a little change. Conventionally, the Fermi-Dirac 
distributions of spin-up and spin-down electron in the electrode are 
different due to finite V s . The higher the temperature is, the less these 
differences are generated. 



In Figure lc and Id, we show I c and I s as a function AT for different 
6 at V g = 0.1 meV and T = 0.02 meV, respectively. In this case, 
interesting phenomena can be observed. When 9 = 0, n, 2%, I c first 
increases and then decreases with decreasing of AT, but it decreases 
monotonically when 9 = n/2. However, J 5 changes monotonically 
with decreasing AT and equals to zero at 9 = 7i/2,3n/2. 

Thermoelectric coefficients. Next, we focus on the thermopower 
phenomena in the linear response regime and assume spin-bias 
only exists at the left electrode. Figure 2a and Figure 2b display 
charge- Seebeck and spin-Seebeck coefficient as a function of 9 for 
different values of V g respectively. At the point V g = 0.1 meV, S c and 
S 5 are zero due to the same weight and the opposite transmission 
direction of the currents and energy carried by electrons and holes. 
At the special points 9 = 0, 2n, the contributions of the spin-up and 
spin -down electrons are the same weight leading to Gf = Gj and S s = 
0, no matter spin bias exists or not. S c and S 5 reach the maximum 
values for each V g when 9 = n. In Figure 2c, we plot all transport 
coefficients as a function of 9 with V g = —0.05 at T = 0.02. One can 
observe that electron thermal conductivity (k e ) is zero at special 9, 
which may cause figure of merit Z C T or Z S T tend to be infinite. 
Furthermore, we plot the k e as function of V g and 9 at T = 0.02 
meV and find that zero k e exists only under special conditions of 
V g and 9 in Figure 2d. It is interesting that one can have ultrahigh 
spin thermopower in a single molecular magnetism through 
manipulating the angle 9 between the easy axis of SMM and the 
spin orientation of the electrodes, and also by tuning V g . 

Thermoelectric coefficients have been investigated through solv- 
ing the non- equilibrium Green's function (NEGF) in detail. The 
general formulas are derived to calculate the currents which depend 
on the angle 9 between the easy axis of SMM and the spin orientation 
in the electrodes and spin bias. The spin bias destroys the SU(2) 
symmetry of electron- spin in nonmagnetic electrodes, which leads 
the fact that the angle 9 influences the redistributions of different 
spin currents. It is amazing that pure spin currents can be obtained by 
tuning V g and AT with a finite spin bias and arbitrary 9 except of 
1 3 

9= -n,-n. In the linear regime, infinite figure of merit can be gen- 
erated by tuning V g at special angles 9 with spin bias. Specially, when 
the angle 9 is equal to zero or 27i, spin thermopowers vanish ident- 
ically even though spin bias exists. These phenomenons may pro- 
vide a new approach for the design of SMM devices based spin 
caloritronics. 

Discussion 

The details of the constituents of the currents at 9 = 0 are shown in 
Figure If. When V g is equal to 0.1 meV, the lowest-energy states of the 
isolated SMM are four-fold degenerate: \5/2 : ±5/2) and 
1 0,2 : +2), and the second-lowest-level state is a double degeneracy: 
1 1,3/2 : +3/2) _ . The energy level difference between the lowest 
and next-lowest levels is 0.0839 meV. The currents are mainly con- 
tributed by the transitions of 1 3/2 : +3/2)"o|0,2 : +2) and 1 1,5/2 : 
+ 5/2) o 1 0,2 : +2). According to Eq. (9), we can approximate 

j ±« ±Y*J1 \ a T 

1 -\\,5/2:±5/2)o\0,2:±2) c,Z ±J\ ^ ) ancl i |l,3/2: ±3/2) _ o|0,2: ±2) °C 

J-L 

/0.0839 + y 5L /2\ /0.0839\ A , , 

/ — I —f — — — I . At the electron-hole symmetry 



Tl J J V T R , 
point, i|i,5/2: + 5/2)o|o,2: + 2) are on ly controlled by V sL and temper- 
ature T L . 

However, I\ lj3 j 2 , ± 3 / 2 ) - <>|o,2: ± 2) depend on the temperatures of the 
two leads, V sL and the energy- difference between the lowest and the 
second-lowest state. From Figure If, one can find that I s is mainly 
contributed by 7| lj5 / 2 . + 5 /2)o|o,2: ± 2) > but I c is decided by all transitions. 
It is amazing that pure spin currents can be obtained at high tem- 
peratures. In Figure lg, it is clear that I c first increases and then 
decreases with the increasing of the average temperature T at V g = 
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T (meV) 



(h) 



T L =0.8 meV, V g =0.27 meV 


V ; =1 OOmeV 




V =0 70 meV 




V_=0.35 meV 




__ V =0 01 meV 




2 0.4 0.6 0.8 1 1.2 1.4 1.6 

T R (meV) 

0.2 meV, / = 0.1 meV, K 2 = 0.04 meV, U = 1.0 meV, 

- T R = 0.0002 meV and T = 0.02 meV: In the (c) and (d), VAs 



Figure 1 | I c and I s as a function of #for different Vg and AT with parameters S = 2, e 0 = 

ic B = 1, P L = P R = 0. In Fig (a) and (b), a tiny temperature difference is considered: AT = T L - 
set to 0.1 me V and the average temperature is fixed: T = 0.02 meV. (e) shows the temperature difference influences on the I c and I s with 0 = 0. 
(f) displays the details of the constituents of the I c and I s with 0 = 0. (g) displays I c and J 5 as a function of Tfor different V^with T L = 0.2 me V and V" 5 = 0.01 



me Vat 0 = 0. (h) shows J c and I s as a function of T R for different V" 5 with T L = 
dash dot lines mark spin currents in (g) and (h). 



0.8 meVand V g = 0.27 meVat 6 = 0. Solid lines denote charge currents and 
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0.27 meV. But I s only has little change and is not equal to zero. Due to 
the spin splitting induced by spin bias, pure spin currents can be 
generated at arbitrary temperatures. In Figure lh, the system tem- 
perature is five times as the anisotropy-induced energy barrier and 
pure spin currents can be obtained by increasing V sL . 

Interestingly, the numerical results show k e can be equal to zero 
with changing V g and 6 in Figure 2d. It means that Z S T(Z C T) may be 
infinite when S C (S S ) and G a (o = t!) are finite. The exact choices of V g 
and 6 are related to the details of the system's parameters. But it is 
necessary that S s must be larger than S c . It is well-known that elec- 
trons move from high temperature to low temperature, and it is not 
related to electronic spin. However, the spin-up and spin-down elec- 
trons under the spin bias can move in the opposite direction, and 
carry different energy according the Eq.10. Due to the competition 
between the temperature difference(AT) and spin bias (V sL ), and the 
scattering between spin-up and spin-down electrons induced by the 
nonlinear spin exchange, it is possible that there are non-zero 
amount for thermopower and electrical conductance when thermal 
conductivity k e is zero. 

Methods 

Non- equilibrium Hubbard Green function has been used to solve the thermoelectric 
transport in the sequential and linear response regime 35 . The system Hamilton can be 
rewritten by the transition operator 36 " 39 , i.e. d a =X 0a + 8 a X d2 with 5 a = + 1( - 1) for 



<7 = TU)(T = j),s + =X^,5_=X^,and5 
operator can be expressed as 40 : 



-X^). For large spin, the 



S+= C+Y m+1 ' m , S~= C m yn ' 

m=-S,S m=-S,S 



S z = Y mY m ^ m . 

m = —S,S 



(2) 



where S is the spin quantum number, Y m,n = \ Sm > < Sn\ and C~ = 
y 7 (S + m + 1)(S+ 1). Finally, the retarded Green Function is written as 

G» = ^ «X 0ff Y fflffl |^ + » r + (5 ff «X s y mffl |^ » r (3) 

By using the Dyson equation and the Keldysh forum, the retarded(advanced) and the 
lesser(greater) Green's function can be compactly expressed as respectively 



_yr/a 
leads, a 



(4) 



(5) 



Here, ^/^^are the electron self-energy in the second-order approximation and 
the formulas for calculation are 



& = + - 

^L / 1,„=+'Et(<+ + <- 

x = L,R L 
d = + - ,k 



(6) 



where 




Figure 2 | Here, we consider spin-bias only exists at the left lead, (a) and (b)show the S c and S s as a function of 9 for different V g respectively. 

(c) displays the thermopowers as a function of 6 with V g = —0.05 meVat T = 0.02 meV. (d) shows conventional thermal conductance (k e ) as a function of 

V g and 6 at T = 0.02 meV. The other parameters are chosen as same as that in Figure 1. 
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A^ = (l + ^P«)rfn 2 (S, 



2 2 



(7) 



(8) 



Here u L } R) = ±(eVj2 + 5 & eV sa /2) is a chemical potential with 



eVj 



l(r) ~- 



and eV< 



2 5L W = - 2 

a' =+( — ). P a denotes the polarization of a lead and is defined as 

Pa = (p a + — P a - )/ 0°a+ + Pa-)- fao' is the Fermi-Dirac distributions of a lead. 

<W' = 1(0) for (7 = |(|) and cr' = + ( — )(o^o'). The voltage is defined as 

V = V L — V^and the spin voltage is written as V s = V sL — y s R.In the Hubbard operator 

representation, the eigenenergies of the unperturbed SMM can be obtained exactly, so 

the equation (5) is rigorous 41 . Following the Landauer-Buttikier formula, the 

expressions for the currents 42 are 



and do, = +l(-l)for 



"' = +,- J 



dcoT a ((D) 



a* + +a*_ a* + +a*_ 



A L a(7 ,fL,c' 



whereT (7 (co) denotes the transmission coefficient of spin-cr electrons 

T a (cn)=r L a G a r a G« 



with 



r" a =r a (i+5„p a cos(e a )). 



(9) 



(10) 



(11) 



(12) 



We can directly obtain the same formula for I* withL<-^K, depicting electric current 25 
l * = if a ~ I 1 and heat currentQ, = (Q L a - Q R a ) /2. I c = J T + Ij and J s = J j — Ij 
denote the charge current and spin current respectively. It is noticeable that our 
formulas (Eq. 9 and 10) are different from the conventional expressions 41 ' 42 , in which 
6 and Fermi-Dirac function are coupled to each other. 

In linear response regime, the thermoelectric coefficients are expressed as 



V + d a Q.5V s 



AT = 0 



V 
AT 

■ E 



= J^( AI i + AI ± 

i c =o,i s =o 2ATUt G i 



i c =oj s= o 



ArVc t 

I S =0J C = 0 



AI, 



(13) 



Here, S c and S s denote the charge Seebeck and spin Seebeck respectively. G a is o- 
electron conductance. k e is the conventional thermal conductance. A/^is thecr- electric 
current induced by a temperature difference at zero voltage bias and zero spin bias. 
AT= T L — T R is defined as temperature difference and the average temperature is 
expressed asT= (T L + T R )/2. Finally, The spin figure of merit 
Z s T = (G j — G i ) S 2 S T/k and charge figure of merit Z c T = (G| + G i ) S 2 C T / k can be 
calculated, k = k e + kp^is the thermal conductance with contributions from both 
electrons (k e ) and phonons(/c^) 43 " 45 . In our model, the phonon transport is not 
considered due to large mismatch of vibrational spectra between the SMM 
and leads. 
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